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Finding genes that cause human hypertension is 
not straightforward, since the determinants of blood 
pressure in primary hypertension are multifactorial . 
One approach to identifying relevant genes is to 
elucidate rare forms of monogenic hypertension. A 
relevant mutation may provide a rational starting 
point from which to analyse the pathophysiology of 
a condition affecting 20% of the world's population. 
In 1973 a family with autosomal dominantly inherit- 
ed brachydactyly and severe hypertension, where 
the two traits cosegregated completely, was 
described 2 . We have now re-examined this kindred, 
and localized the hypertension and brachydactyly 
locus to chromosome 12p in a region defined by 
markers D12S364 and D12S87. As the renin- 
angiotensin-system and sympathetic nervous sys- 
tem respond normally in this form of hypertension, 
the condition resembles essential hypertension. 
This feature distinguishes this form of 
hypertension from glucocorticoid 
remediable aldosteronism and 



ger was shortened. When a ratio of index to middle fin- 
ger length was plotted, affected and nonaffected persons 
were separated completely (Fig. la). Therefore, brachy- 
dactyly was selected as the phenotype for linkage analy- 
sis. The relationship between median arterial blood 
pressure and age for the two groups shows a marked dif- 
ference in slope (Fig. 2b). There was an observed 50 mm 
Hg difference in blood pressure between the two groups 
in the fifth decade of life. Historically, the cause of death 
in affected individuals was stroke below age 50 years. We 
excluded children <10 years in this particular analysis, 
since their phenotypes were not invariably discernible. 
Affected and nonaffected persons , did not differ with 
respect to body mass index. Their years of schooling 
and literacy status were not different; therefore, affected 
persons were not mentally retarded. 

We first performed an exclusion analysis in a subset 
of 36 family members with -300 DNA markers. Only 
markers on chromosome 12 suggested linkage to the 
brachydactyly phenotype. A maximum pairwise lod 
score of 2.62 was obtained with marker D12S83 at a 
recombination fraction 0 = 0.10. The complete family 
(n= 55) was then used to test for linkage with addi- 
tional markers from this region (Table 1). The critical 
chromosomal region was defined by the two polymor- 
phic markers D12S364 and D12S87 which are approxi- 
mately 24 cM apart and reveal obligatory 
recombination events. In multipoint linkage analysis, 
the lod score reached a peak value of 9.29 at the posi- 
tion of marker D12S310 (Fig. 3). 

Two autosomal dominant forms of hypertension 
have been identified so far. In glucocorticoid remedia- 
ble aldosteronism a chimaeric gene incorporating fea- 
tures of the 11 G-hydroxylase and the aldosterone 
synthase gene was identified 3 " 5 . In Liddle's syndrome 
mutations in the and y subunits of the epithelial sodi- 
um channel were demonstrated 6 ' 7 . Elucidation of these 
conditions has contributed greatly to under- 
standing basic mechanisms and has 
also enhanced patient care 1 . 
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Uddles syndrome, which 
are salt-sensitive forms of 
monogenic hypertension 
with very low plasma 
renin activity 3 " 7 . We 
suggest that identifi- 
cation of the gene 
involved in hyper- 
tension and brachy- 
dactyly and its 
mutation will be of 
great relevance in 
elucidating new 
mechanisms leading 
to blood pressure ele- 
vation. 

In this family hyperten- 
sion and brachydatyly show an 
autosomal dominant mode of 
inheritance (Fig. 1). All affected per- 
sons had brachydactyly 
mainly involving the 




To our knowledge, our fami- 
ly 2 provides the third exam- 
ple of autosomal dominant 
hypertension. Pheochro- 
mocytoma and inherit- 
ed renal disease have 
been excluded 2 . The 
short, fourth 
metacarpal of some 
subjects resembles the 
anomaly associated 
with Albright's heredi- 
tary osteodystrophy. 
However, our subjects 
were neither obese nor 
mentally retarded. Their 
calcium and phosphate 
homeostasis was normal and we 
found no linkage on chromosome 
20 (ref. 8). An Albright's osteodystro- 
phy variant with 

Fiq 1 Hypertension and brachydactyly pedigree. Twenty-five affected {24 brachydactyly has 



" • livina and 1 since deceased,' Tv.2) 7 and 21 unaffected members were reC ently been mapped 
metacarpals to varying " v "' a ^ ™« 



degrees (Fig. 2a). 
Notably, the index fin- 



nvtnq ana l since ueiw<u>eu, 'v.«-/ 

included in the linkage analysis. In addition, 9 children with questionable 
Dhenotype status (?) were included. Affected individuals are represented 
as hatched figures. A diagonal line through symbol indicates decerased. 



to chromosome 2q37; 
however, these indi- 
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Fig. 2 a, Brachydactyly (typical example 
left) as demonstrated by the quotient 
between the index finger and the middle 
finger, b, Median arterial blood pressure 
differences between affected (solid figures) 
and nonaffected (open figures) family 
members. The two individuals marked in 
brackets received antihypertensive treat- 
ment. 
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viduals were retarded and not hypertensive 9 . 

An additional feature of affected persons in the 
hypertension brachydactyly kindred was an increase in 
the rate of fibroblast growth compared to age and gen- 
der-matched controls (unpublished observations). 
Spontaneously hypertensive rats also exhibit an 
increased rate of fibroblast and smooth muscle cell 
growth compared to Wistar-Kyoto control rats 10 . This 



feature underscores the resemblance of this inherited, 
form of severe hypertension with the closest genetic 
model to essential hypertension in man. 

We believe that the two phenotypes, hypertension 
and brachydactyly, are caused by either a single 
pleiotropic gene or two closely situated genes. Cytoge- 
netic analysis showed no aberrations at 850 band reso- 
lution. Potentially interesting candidate genes include 



Table 1 Pairwise lod scores between brachydactyly and hypertension and chromsome 12 marker loci 



Locus 


Alleles 


PIC 


0.00 


D12S336 


5 


0.65 


— oo 


D12S320 


8 


0.78 


— oo 


D12S364 


11 


0.81 


— oo 


D12S308 


3 


0.51 


! 1.671 


D12S310 


6 


0.56 


11.806 


D12S363 


5 


0.68 


5.247 


D12S799 


11 


0.81 


11.439 


D12S87 


8 


0.69 


— oo 


D12S345 


9 


0.81 


—oo 


D12S59 


11 


0.80 


—co 


D12S85 


6 


0.72 


— oo 


D12S361 


6 


0.69 


— oo 



Recombination fraction (0) 






©max 


0.01 


0.05 


0.10 


0.20 


0.30 


^max 


-1.781 


0.680 


1.435 


1.644 


1.237 


1.672 


0.170 


3.138 


4.134 


4.189 


3.499 


2.350 


4.221 


0.079 


1.276 


2.501 


2.753 


2.433 


1.686 


2.755 


0.106 


1.637 


1.499 


1.325 


0.967 


0.603 


1.671 


0.000 


11.309 


10.553 


9.569 


7.456 


5.031 


11.806 


0.000 


5.188 


4.906 


4.483 


3.486 


2.360 


5.247 


0.000 


11.250 


10.479 


9.475 


7.328 


4.979 


11.439 


0.000 


4.068 


4.890 


4.750 


3.802 


2.519 


4.901 


0.059 


-0.076 


2.253 


2.848 


2.748 


2.033 


2.925 


0.135 


4.067 


5.515 


5.574 


4.667 


3.200 


5.626 


0.077 


2.626 


4.142 


4.299 


3.636 


2.488 


4.313 


0.087 


3.827 


4.706 


4.630 


3.777 


2.556 


4.733 


0.065 



cM IV.24 UL6 V.42 V39 HL3 V.15 VJ 




Ra 3 a Multipoint linkage map using the most informative markers, b, Schematic presentation of chromosome 
n IriS the location of the hypertension brachydactyly locus with defined flanking makers. Sagging 
haplotype of subject IV.24 and haplotypes of other subjects showing am « man ™*£*™> ^ 
area represent the minimal cosegregating region. Numbers correspond to those in the family tree (Rg. 1). Solid 
figures represent affected persons, while the open figure represents a nonaffected individual. 
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the gene for parathyroid hormone related peptide 
(PTHrP). PTHrP is expressed by vascular cells, lowers 
blood pressure, and attenuates cell proliferation 11 . A 
recent gene disruption experiment resulted in 
homozygous mice with skeletal anomalies 12 . However, 
we have since ruled out the PTHrP gene as the gene 
responsible for this syndrome and are continuing our 
search. In our subjects a single gene is responsible for a 
50 mm Hg elevation in blood pressure by age 50 years. 
Elucidation of this gene will likely give new insight 
into mechanisms of blood pressure elevation and may 
be of relevance to primary hypertension. 

Methods 

Family material. We visited the family described by Bilgintu- 
ran et al}, which lives in a remote area on the north-eastern 
Black Sea coast of Turkey and examined 45 members after 
written, informed consent was obtained. Six additional mem- 
bers currently residing in Germany were examined 4 weeks 
later. Nine additional members were identified in a second visit 
to Turkey. Height and weight were obtained, as well as pho- 
tographs of each individual's hands. Blood pressure and heart 
rate were measured by an oscillometric automated method. 
Venous blood was obtained for DNA extraction. 

DNA-marker analysis. Genotyping was performed using the ABI 
PRISM Genotyping System, including the Linkage Mapping Set, 
PCR 9600 thermocyclers, ABI DNA Sequencers, GeneScan and 
Genotyper software from Applied Biosystems Division, Perkin- 
Elmer Corporation (ABD). The PCR primers contained in the 
Linkage Mapping set amplify dinucleotide repeat loci spaced at 
approximately 10 cM. The loci were selected from the human 
linkage map generated by G6nethon 13 based on map position, 
heterozygosity and performance in routine PCR analysis. The 
markers in the Linkage Mapping Set are organized in 28 panels, 
with 9 to 17 primer pairs per panel whose products can be elec- 
trophoresed and detected in a single gel lane. Forward primers are 
labelled with either 6-FAM, HEX, or TET fluorescent dyes which 
can be distinguished due to their different spectral properties. 
PCR products related to specific markers are automatically identi- 
fied in each lane by their molecular weight and fluorescent dye 14 . 
PCR reactions and electrophoresis were performed according to 
manufacturer's protocol. Electrophoresis and detection were done 
on both, an ABI 373 DNA Sequencer equipped with GeneScan 1.2 
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software and an ABI 377 DNA Sequencer equipped with GeneS- 
can 2.0 software. Data were exported as a text file from Genotyper 
for subsequent linkage analysis. 

Linkage analysis. Prior to linkage analysis, we performed exten- 
sive simulation studies using the SLINK simulation program 15 . 
Expected lod scores were determined for various subsamples of 
the pedigree under different genetic models for both, the trait 
and the marker locus, and at recombination fractions ranging 
from 0.01 to 0.1. The pedigree subsample chosen for the genome 
scan was 36 family members (21 affected, 15 nonaffected), and 
yielded an E-lod of +4.85 for a codominant marker locus having 
a heterozygosity of 0.75 and linked to the fully penetrant trait 
locus at a recombination rate of 0.05. Linkage analysis was per- 
formed using the LINKAGE package version 5.1 (ref. 16). For the 
trait locus we used a model of nearly complete penetrance (0.99) 
and low phenocopy rate (0.001), a disease allele frequency of 
0.001 and absence of mutation. Children below age 10 years were 
typed as unknown. The codominant marker locus was individu- 
ally modelled according to allele numbers and frequencies. Once 
significant linkage to chromosomal 12 was found, we switched to 
the full pedigree, including one loop (person TV. 14, see Fig.l), 
and 9 children typed unknown for the trait phenotype and thus 
providing marker information only. This full system was used to 
calculate two-point and multi-point lod scores, using the 
FASTLINK package 17 . A multipoint map was constructed with 
the pedigree divided into three parts and with the four most rele- 
vant marker loci (allele numbers down-coded) using LINKMAP. 
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In 1992, we described a second-generation genetic linkage map of the human genome. j 

Using 1,267 new microsatellite markers, we now present a new genetic linkage map | 
containing a total of 2,066 (AC) n short tandem repeats, 60% of which show a 
heterozygosity of over 0.7. Statistical linkage analysis based on the genotyping of eight 

large CEPH families placed these markers in the 23 linkage groups. The map includes j 

1 266 intervals and spans a total distance of 3690 centiMorgans (cM). A total of 1 ,041 j 

markers could be ordered with odds ratios greater than 1 000:1 . About 56% of this map j 

is at a distance of 1 cM or less from one of its markers. j 
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The use of microsatellite markers 1,2 has permitted the 
construction of new genetic maps of the human genome 
as well as those of other mammals 3 " 5 . These maps can be 
used to map any Mendelian trait and in particular, 
monogenichuman diseases. On theother hand, theaverage 
length of the intervals between adjacent markers often 
requires a great deal of work to isolate new markers which 
are even closer to mapped genes. In order to facilitate this 
work for the entire human genome and to refine the 
genetic linkage intervals that can be covered rapidly by 
cloned DNA fragments for gene identification, we have 
developed and mapped 1,267 new (AC) microsatellite 
markers. These new markers have been integrated into 
the map which we constructed previously and described 
at the end of 1992 (ref. 4). 

New markers 

New markers were obtained using the same procedures as 
for our first 814 markers 4 . However, as there were some 
indications that subtelomeric regions were less densely 
covered with markers 4 , 62 of the markers on this new map 
(markers AFMal20to AFMal52) are from the H3 and H2 
isochores which are preferentially located in terminal 
parts of chromosomes 6 . Although the chromosomal 
distribution of the current complete collection of markers 
isverysimilartothatobtainedpreviousl/^hedistribution 
of markers from the H2 and H3 isochores is quite different 
(Wunderle et al y manuscript in preparation). 

The distribution of heterozygosity of the 2,066 markers 
is shown in Fig. 1. Compared to the previous map of 814 
markers, which had a mean heterozyosity of 0.75, the 
average heterozygosity is 0.70, and the percentage of 



markers with a heterozygosity of over 0.70 has decreased 
from 74.4% to 60.1%. This reduction is due to the use of 
all the markers which showed at least three alleles when 
tested on four indivi duals 4 , whereas in the 814-marker 
mapwe selected essentially those markers for which the 
frequency of the majority allele was less than 0.5. 

Genotyping and map construction « 

Genotyping was carried out according t©4h£-m*ikipi£2L 
procedure previously described 7 . Although we tried to 




-4" 

0.2 0.3 0.4 0,5 0.6 0.7 

Heterozygosity 

Fig. 1 Comparison of relative distribution of heterozygosity 
values of the (CA) n microsatellites in the 1992 and 1993-94 
versions of the Genethon map. Dark columns: 814 markers; 
H > 0.5 for 97% and H > 0.7 for 74.4%; mean 
heterozygosity value = 0.75. Open columns: 2,066 markers; 
H > 0.5 for 93% and H > 0.7 for 60.1 %; mean 
heterozygosity value = 0.70. 
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avoid genbtyping microsatellites which were identical to 
sequences deposited in databases by others, occasionally 
the sequences of some of our markers were deposited in 
databases while our work was in progress. We have retained 
these markers on our maps because they are useful as 
common reference points among several maps (see notes, 
Fig. 2). 

Genotyping errors were detected using the same 
diagnostic software programmes as described previously 4 . 
Genotypes of markers which appeared as double 
recombinants after several verifications were used for 
construction of the maps, although most of these double 
recombination events probably correspond to mutations 
or gene conversions. 

A total of 305 mutations was observed in 278,338 
genotypings. This percentage of 0.1% (0.05% of haploid 
genotypes) is very close to that observed previously by us 4 
and by others 8 . About half of these mutations probably 
result from the use of DN A from continuous lymphoblastic 
cell lines 8 . 

The procedure used to construct the map was again 
similar to the one used for the previous version. However, 
our published map of 8 14 markers was used as a framework 
to position the new markers. Several markers which 
appeared on our previous map have been removed. In the 
majority of cases (AFM263zh9, AFM046xcl 1, 
AFM224xflO, AFM234tg3, AFM200ya9, AFM163xal, 
AFM238ybl0 and AFM182xg9), these markers had one 
or more alleles which could not be amplified with the 
primer pairs that we selected (cryptic alleles). These 
markers will be returned to the map when the primers 
have been modified so that all the alleles can be amplified. 



• AFM136xe3, AFM165zd4, AFM217yc5 and AFM120xf6 
each corresponded to two polymorphic loci and were 
therefore eliminated. During the initial sequence 
comparison, it was not noticed that marker AFM186xd2 
was identical to AFM070ya9, so that the former was 
elimininated. Marker AFM262vg9 had a tendency to 
produce nonspecific bands and was also removed from 
the map. Finally, marker AFM123yf8 could not be 
repositioned unambiguously on the present map. 

The new map 

The new map covers a total distance of 3,690 cM, which 
represents an increase of 114 cM over the 814-marker 
map (Table 1 ). The extensions, which account for 135 cM, 
are due to the addition of more telomeric markers (Table 
1 ) . The total increase is less than these telomeriqextensions 
because of modifications of distances in the interior regions 
of some chromosomes. These modifications usually 
involve small decreases in distances, but also increases. 
No new markers were added to the ends of eleven 
chromosomes: 3-10, 15, 18 and 19. However, chromosome 
6 and especially chromosome 3 show significant increases 
in their total length, whereas chromosome 4 shows a 
marked decrease. Chromosomes 1 and X show decreases 
in total length despite the addition of one or more markers 
to at least one of their extremities, and 14 and 17 show 
increases which are clearly less than the extensions of their 
extremities. Sex-specific distances of each chromosome 
are indicated in Table 1. 

There are a total of 1,266 intervals (Fig. 2). in the map 
which corresponds to an average distance of 2.9 cM 
between markers. We were able to position 1 ,04 1 markers 



Table 1 Comparison of the main features of the 1992 (ret 4) and 1993-1994 versions of the Genethon hum an genetic linkage map 

Chromosome Length covered (cM) 



1 



y 

4 
's; 



l 3 • 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 

1 13q 
s 14q 
h 15q 
i? 16 
? 17 
k 18 
l f| 19 
t 20 ' 

21 q 
I 22q 
% X 

i'C. Total 



1992 
sex- 
average 

295 
277 
221 
229 
201 
201 
195 
155 
160 
178 
161 
172 

99 
125 
107 
119 
128 
126 

95 
101 

29 

32 
170 

3576 



1993-94 
sex- male female 
average 



Largest 
gaps (cM) 

1992 1993-94 



Map extension Number of Number of markers 

since 1 992 (cM) markers mapped positioned with 

odds > 1000:1 



1992 1993-94 



1992 1993-94 



pter or 
qcen 



qter 



292 
281 
235 
211 
201 
208 
196 
155 
158 
179 
160 
180 
123 
127 
107 
130 
133 
129 
99 
120 
49 
49 
168 



218 
211 
193 
150 
151 
130 
131 
99 
121 
136 
119 
138 
102 
103 
80 
101 
108 
99 
86 
83 
45 
40 



362 
334 
282 
263 
251 
274 
246 
202 
189 
212 
183 
212 
146 
153 
135 
164 
155 
161 
119 
142 
61 
67 
168 



19 
24 
■ 19 
17 
14 
20 
'15 
18 
18 
13 
14 
30+37 
12 
17 
14 
27 
15 
25 
32 
18 

9 

8 

22 



14 
11 

8 
10 

7 
13 

8 
18 
14 
12 

8 
17 
19 

9 
12 
13 
13 
13 
24 
10 
12 

8 

18 
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with an odds ratio of 1,000:1 or better. These maps, 
constructed with several tools based on the LINKAGE 
programme package, were compared to those using the 
same genotypes processed with the MultiMap algorithm 9 , 
based on CRI-MAP 10 . Both sets of maps were essentially 
identical in order as well as distance. The main difference 
was observed with a few markers that were rejected by the 
GMS algorithm (Gene Mapping System) 11 because they 
could not be mapped to a region that was sufficiently 
precise whereas they were included in the comprehensive 
maps resulting from the MultiMap process. Refinement 
of tie map and the increase in the number of markers has 
led to some modifications in the previous order of the 
markers. Four of these modifications concern markers 
previously ordered with odds ratios of greater than 1 ,000: 1 
and 14 concern markers positioned with lower 
probabilities. 

Distribution of markers 

Only one gap of over 20 cM remains on the map. The 
other gaps have been reduced in size. There are only 22 
remaining gaps of over 10 cM, which represent 6 of the 
gaps which were over 20 cM on the previous map. A 
significant proportion of the markers from the H2 and H3 
isochores (8 out of 62) were found to map to the distal end 
of the chromosomes and a number of the others are 
subtelomeric, as expected 12 . This indicates that markers 
from these GC-rich regions should permit a more dense 
coverage of numerous subtelomeric regions. 

Correspondance between genetic and physical distance 
must await integration of genetic and physical maps. 
Genetic linkage maps integrating polymorphic markers 
from different sources including AFM markers from the 
first set of 814 have been established recently 9,13 . A more 
extended integration project using a different strategy 
and including new markers from the present map is in 
progress. 

Conclusion 

About 56% of our latest genetic linkage map of 3,690 cM 
is at a distance of 1 cM or less from one of its markers. 
In many cases, these distances can be covered by cloned 
DNA sequences 14 . Moreover, the isolation and mapping 
of 3,000 additional markers is in progress. This will increase 
the density of marker coverage and perhaps extend some 
of the chromosomal maps. This should accelerate 
considerably positional cloning of hereditary disease genes 
by facilitating the search for additional close genetic 
markers and candidate exons. 

Methodology 

Marker development. Marker development was carried out 
essentially as described*. DNA libraries were made from an A/uI 
DNA digest of 46>XX human DNA (sized between 300-500 bp) and 
cloned in Ml 3. The sequences of the templates from the (CA) n or 
(GT) positive clones were used to define PCR primers. The 
synthesized primers were tested on four unrelated 46,XX individuals 
to obtain a first estimate of the polymorphism of the tested 
microsatellite markers. Markers with three or more alleles were first 
assigned to their chromosome and genotyped as described 4 . 
The H2 and H3 isochores were isolated by caesium sulphate 



density gradient centifugation of total human DNA in the presence 
of BAMD 6 . The (AC) n markers found in the fractions with the 
highest GC content were isolated by a single Alul digestion shotgun 
procedure as described 4 . 

Genotyping. Individuals from the eight CEPH families (102, 884, 
1331, 1332, 1347, 1362, 1413 and 1416) were genotyped using 
standard procedures, as described'. PCR amplifications were 
performed in 50 ul reaction mixtures, containing 40 ng of genomic 
DNA, 50 pmol of each primer, 125 \M dNTPs, 10 mM Tris pH 9, 50 
mMKCl, 1.5 mM MgCl 2 , 0.1% Triton X- 100, 0.01% gelatin and 1 U 
of Taq polymerase (Amersham). Amplifications were carried out 
using the "hot-start" procedure, in which the Taq polymerase is 
added to the reaction mixtures after a first denaturation step (5 min 
at 96 °C) after which 35 cycles of denaturation (94 °C for 40 s) and 
annealing (55 °C for 30 s) are performed. An elongation step (2 min 
at 72 °C) ends the process. For each DNA sample, 16 amplification 
products from different markers were ethanol- precipitated and 
loaded together into single lanes of 6% polyacrylamide-8M urea 
denaturing gels. After migration, the DNA was transferred from the 
gel to a Hybond N + nylon membrane (Amersham) by a contact 
blotting procedure 7 . 'The markers were then revealed by successive 
hybridizations with one of the PCR primers which was peroxidase- 
labelled by modification of the ECL procedure (Amersham) and 
exposed to autoradiographic X-ray films 7 . ► 

Map construction. Markers were assigned to chromosomes by 
pairwise linkage and possible genotyping errors were identified by 
comparisons between families of the pairwise recombination events 
between linked markers. After genotype corrections, markers from 
a chromosome-specific dataset were positioned on a framework 
consisting of the map of 814 markers 4 using a map construction 
algorithm. The order of markers in the framework and complete 
maps were determined with the GMS algorithm 11 . Briefly, 
recombination estimates for a preliminary, or trial order of the loci 
are used to divide the loci into subgroups of closely linked loci. 
Likelihoods are evaluated for different placements of subgroups, and 
for alternative' orders of the loci within each subgroup. The best- 
supported order (i.e. the order with the greatest likelihood) is chosen 
as new trial order, and iterations are continued until convergence. 
Based on the best -supported order for the framework map, 
recombination fractions between adjacent markers were estimated 
with the LINKAGE programs 15 . Markers from this framework that 
underwent corrections since the 8 14-marker map were processed as 
new incoming markers. This led to a provisional order which was 
further reassessed as described 16 . Once the order remained unmodified 
after further computation, a search for double recombination events 
was undertaken. The maps were re-evaluated using the corrected 
genotypes until no further double recombination event could be 
eliminated. 
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